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Abstract. 

The interplay between chaotic tunneling and dynamical localization in mixed phase 
space is investigated. Semiclassical analysis using complex classical orbits reveals that 
^SJ I tunneling through torus regions and transport in chaotic regions are not independent 

processes, rather they are strongly correlated and described by complex orbits with 
both properties. This predicts a phenomenon analogous to the quantum suppression of 
classical diffusion: chaotic tunneling is suppressed as a result of dynamical localization 
in chaotic regions. This hypothesis is confirmed by numerical experiments where the 
effect of destructive interference is attenuated. 
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Tunneling phenomena are purely quantum effects. Nevertheless its nature is 
strongly inffuenced by underlying classical dynamics. In particular nontrivial aspects 



CN ! appear in the dynamical tunneling in mixed-type phase space, in which quasi-periodic 

Q ! and chaotic trajectories coexist [1]. Tunneling transitions between quasi-doublet states 

! are enhanced by chaotic states [2] and the existence of nonlinear resonances also 

^ ! leads to a qualitative change of tunneling processes [3,4]. The problem of quantum 

! tunneling in multidimensional systems or more specifically in nonintegrable systems 

^ ' first raised in [5] would be of fundamental importance in quantum mechanics, and an 



approach taken there has recently been extended in [6]. However, our understanding for 
multidimensional tunneling is still far from complete. A primary difficulty lies in the 
5^ ! fact that dynamical tunneling in chaotic systems takes place in very complicated phase 
space; the structure of classical phase space itself is not an easily understandable object. 
Another reason would be that dynamical tunneling proceeds in complex environments 
in the sense of wave phenomena. Scarring [7] or dynamical localization [8,9], or other 
types of invariant structure become sources of partial structures often observed in wave 
functions. It is not clear at all to what extent these various wave effects are independent 
of each other. These aspects make it difficult to evaluate the tunneling rate between 
torus to chaotic regions quantitatively. 

The trajectory description would be one of promising strategies to understand 
quantum phenomena of chaotic systems. In particular, the semiclassical analysis is now 
recognized as an efficient approach to this end, and there are indeed a bunch of numerical 
tests supporting its validity. Concerning dynamical tunneling, it was shown that the 
complex semiclassical theory works fairly well and explains the mechanism of tunneling 
penetration out of quasi-periodic regions to chaotic seas [10-12]. In such a treatment, 
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the classical dynamics is extended to complex phase space, and the trajectories on the 
Julia set is most responsible for reproducing tunneling wavefunctions [13]. 

The aim of the present Letter is, based on the arguments predicted by the complex 
trajectory description of chaotic tunneling, to show that dynamical tunneling in mixed 
phase space and dynamical localization are strongly correlated to each other, so that the 
destruction of coherence, or more precisely the destruction of destructive interference 
in chaotic regions, not only induces delocalization of wavefunction in chaotic regions, 
but also causes strong enhancement of tunneling transition. In other words, "genuine" 
chaotic tunneling is suppressed by dynamical localization in surrounding chaotic regions, 
which is entirely analogous mechanism as dynamical localization suppresses classical 
diffusion [8,9]. The result must serve as further understandings of amphibious states 
[14, 15], recently discovered quantum states that show the failure of semiclassical 
wavefunction hypothesis [16]. 

We first provide evidence for why we can predict that chaotic tunneling is 
tightly correlated with the dynamical localization process. As found in [11, 13], an 
exponentially large number of complex orbits appear as the contributors of the time- 
domain semiclassical propagator, and they indeed have almost equal weights in the 
semiclassical sum, which means that chaotic tunneling occurs as a consequence of 
superposition of exponentially many component waves. To demonstrate it, we here 
employ an area-preserving map 

where 

T'{p) = ap+ -{di- d2) + - [ap-u + di] tanh6(p - pd) 

+ ^ [-ciP + u; + 6/2] tanh h{p + pd) (2) 

V\q) = -Ksinq. (3) 

As illustrated in Fig{T]^d), phase space is divided into quasi-periodic and chaotic regions 
for 6 3> 1. The kinetic term is almost linear for |j9| < pd, and tends to the standard 
map for \p\ > pd- The parameters di and c?2 was put in order to get rid of small island 
structures which may appear around the border between torus and chaotic regions. The 
smoothing factor tanh is introduced to allow analytical continuation of the map into 
complex plane. 

The wavepacket launched at p = goes out of the torus region due to tunneling 
effects. Time evolution is described by the propagator in p-representation {pn\U\po) and 
its leading-order semiclassical approximation takes the form as 

U'\pQ]Pn) = ^v4fc(po,Pn)exp |^5'fc(po,Pn)| , (4) 

where the summation is taken over all classical paths k satisfying given initial and final 
momenta, Ak{po,Pn) and Sk{po,Pn) stand for the amplitude factor associated with the 
stability and the corresponding classical action, respectively. 
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For given a G M, we have a set for initial conditions of semiclassically contributing 
complex paths as 

= {(Po, go = e + e I po = a, Pn e M }. (5) 

As shown in Figs{T](a) and (b), there exist a bunch of complex trajectories that satisfy 
the boundary conditions imposed on the initial and final steps. The set is composed 
of strings, each of which correspond to the contribution k in the sum (jlj). In case of the 
present map, the number of strings is infinite even for a finite fixed time step, reflecting 
that the map contains transcendental functions, while only finitely many complex orbits 
appear in case of the map composed of polynomial functions. What was found in [11,13] 
is, irrespective of the form of maps, that the tunneling wavefunction is well controlled 
and quantitatively reproduced by special types of complex orbits. Such orbits appear 
as a chained structure in the set A^°, an example of which is displayed in Fig{Tt^b). 
The orbits forming chained structures are shown in red color. It was also shown 
in [11, 13] that the number of complex orbits substantially controlling the tunneling 
process increases as a function of time. Further important features actually shown here 
is that all these approach the real plane exponentially. These are both explained by the 
fact that such orbits are attracted by unstable periodic orbits on real plane [24]. All 
the stable manifolds of unstable periodic orbits extended to complex plane in general 
reach and intersect an initial state {po,qo = ^ + irj) G C^, thereby the orbits launched 
at the initial state go out of regular regions by following stable manifolds and approach 
the real plane exponentially. 

It is important to note that, after reaching quasi- real regions, the dynamics of 
complex orbits is almost governed by the real dynamics. Figures [U^c) and (d) plot the 
distance from the real plane, and the projection of final manifolds onto the real plane, 
respectively. Three curves shown in FigHJ^c) show the itinerary of complex orbits staring 
at three initial points marked in Fig{T](a). Note that each one is not the itinerary of a 
single orbit, but a representative of exponentially large number of orbits around each 
initial point, since these initial points are very close to each other. For example, the 
orbits starting at the red part in Fig{T](b) behaves in the same way as the red curve in 
FiglD^c), so these are not distinguishable from each other. 

Along the final manifold presented in Fig{T](d), |Img„| is very small, meaning that 
it almost follows the stretching and folding mechanism in real phase space. We notice 
that, within a single step, the orbits already go out from the torus region, and then go 
down to the real plane exponentially with being stretched in the unstable direction. In 
the final step (ra = 5 in this case), the stretched manifold is almost real. Therefore, if 
one focuses on the manifold contained in the chaotic regions, the situation is almost the 
same as what is taking place in real dynamics. In this way, complex orbits controlling 
chaotic tunneling bear an amphibious character: running through the torus region in 
purely imaginary space, and extending over the quasi-real chaotic region. 

Now, recall that the wave function in chaotic regions is dynamically localized [8]. 
Although the semiclassical interpretation of dynamical localization is not still clear [17], 
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Figure 1. (a) A set of initial conditions Ai^ (o^ = 0)i a-nd (b) its magnification. Red 
lines are parts of chained structures which substantially contribute to the semiclassical 
propagator |Inig„| of such orbits are very small. The parameters are chosen as 
n = 5, a = 5, 6 = f 00, di — —24, d2 — —26, uj — l,pd = 5 and K = 2. (c) The distance 
from the real plane (D„ — y^\pn\^ + |9np) as a function of time step, where Pn and (?„ 
denote (p, g) at n. The initial conditions are taken from the points indicated in (a). 
Note that the orbits whose initial imaginary parts are small (yellow and green) are not 
necessarily go to the real plane directly, rather take a side trip in complex phase space, 
(d) The Lagrangian manifold projected onto real phase space. Here the Lagrangian 
manifold refers to F"(Af "), and the initial conditions giving these manifolds shown in 
(d) are red parts shown in (b). The green dot (labeled as n = 0) represents the initial 
condition and black curves (labeled as n = 5) is the final Lagrangian manifold. 



it may be at least true that, if the semiclassical description works, the orbits contributing 
to the semiclassical sum (jlj) should have appropriate correlations among them [18], 
otherwise, classical diffusion will be restored since the random phase in the semiclassical 
sum @ cancels the off-diagonal contributions. In the present situation it would be 
natural to expect that the quasi-real orbits discussed above are correlated in chaotic 
regions as well. 

To check this speculation, we add noise in the chaotic region and destroy the 
coherence predicted here. The experiment is done by putting the noise term to the 
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kinetic term Tq{p) as 



T'(p) = r^(p)+eTV,,(p). 



(6) 



Here e represents a stochastic variable obeying the Gaussian distribution with a variance 
cr. As schematically shown in Figj2](a), setting T^oisc,{p) = 1 for |p| > L, = for |p| < L, 
we apply noise only in the region \p\ > L. To see the tunnehng amplitude, we will 
monitor the probability P^°^^ defined as 



where ipn{p) represents the wavefunction at time step n, which is initially placed at 
Po = 0, and Pa and pb are the coordinates specifying the minimum and maximal values 
of the torus region projected onto the p-axis (see also FiglT](d)). Here, we imposed the 
periodic boundary condition in the p-direction. As shown in Figl2]^a), wavefunction 
outside the torus region spreads as time proceeds and the profile of its tail shows that 
the classical diffusion process is recovered due to the destruction of interference which 
causes dynamical localization [19,20]. Figure [2]^b) plots -P^°™**, clearly demonstrating a 
drastic enhancement of tunneling probability. Note that the ensemble of real classical 
orbits whose initial distribution is set to be the same as the quantum initial distribution 
almost stays inside the torus region and does not leak out even under the same noise. We 
can see that with the increase of L, which makes the chaotic region surrounding the torus 
region large, the tunneling rate decreases. As will be reported elsewhere [23], extensive 
computations reveal that the exponent obtained by fitting an exponential function 
to initial decaying interval of -P*°™^, up to 10'' or 10^ in case of Figl2]^b), decreases 
monotonically as a function of L. These are consistent with our interpretation that the 
presence of localized regions suppresses potentially existing tunneling amplitude, which 
is related to an exponentially large number of complex orbits shown in FigUl At the 
same time, note that the noise average is not neccesarry for the enhancement of the 
tunneling: Fig. 2 (a) and (b) are such examples. 

An important remark would be that the tunneling leakage continues but its rate 
7 could not be fitted by a simple exponential function in the whole time scale. The 
rate of penetration slows down with time and it takes very long time, if so, to 
reach the equilibrium state. If the wavepacket spreads over phase space equally, the 
final probability should take the value Pgq'^^ = 0.016, which is estimated under the 
assumption that the wavepacket is uniformly distributed over the phase space. (Note 
again that phase space in the experiment is compact since periodic boundary conditions 
are imposed on both directions.) However, P^°''"'^ is still far above P^q™'^, and the 
localized peak around the torus region remains, meaning that the wavepacket does not 
spread over phase space in an equal weight even after sufficiently long time. 

Since an exponentially large number of complex paths exist also from the chaotic 
region to the torus region [13], the inverse tunneling process should be enhanced as 
well. Indeed, as shown in Figj2](c), if we place the initial wavepacket in the chaotic 
sea, instead of the torus region, the tunneling flow into the torus region is also highly 
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Figure 2. (a) Time evolution of wavepacket launched at po = state. Noise 
is applied only in the hatched blue regions in the schematic picture of phase 
space, a = 4 + 0.1\/3, u = 27r(l + \/7)/10, di = —apd + uj, d2 = —apd — u), 
Pd = TT, L = 27r, \pa — Pb\ = l.Ovr, N = 1000 and h = O.lvr, where N denotes 
the system size. The noise intensity is given as e = O.lvr. (b) P^°^^'^ as a 
function of time step n. The parameter values are the same as (a). Note 
that Peq^^ = 0.016 is far below a saturated value for the case L = 2tt. (c) 
The probability inside the torus region as a function of n for the cases 

without noise (e = 0.0) and with noise (e = 27r/100). The initial wavepacket 
is set as tpo{p) = Aexp{27rir]{p)) (for \p\ > 27r), = (otherwise). Here A is 
a normalization constant, and r](p) = uniformly random in [0, 1]. The plot is 
given after averaging over 10 ensembles with respect to the phase ??(p). 



enhanced with noise being applied on chaotic seas. This result tells us that the chaotic 
region with external noise does not work merely as a sink or reservoir. If this is the 
case, the unidirectional flow from the torus to chaotic regions is expected to take place 
and the tunneling back process should not be observed. 

Similar strong enhancement is observed when we design the system so that the 
wavepacket moves freely in the regions \p\ > L. The motivation to put ballistic regions 
is again in order to get rid of dynamical localization, as the case where noise is applied. 
To realize it, we replace the kinetic term ([2]) by 

T'{p) = -ap{ta.nh.b{p + L) — tanh6(p — L)} 
+ ]-{ap - cj){tanh6(p - Pd) - tanh b{p + pd)}. (8) 
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Figure 3. (a) The wavefunction |^n(p)P at a relatively short (n = 3000:red) 
and a sufficiently long time step (n = 2 x 10^:green) in case of the map with 
the kinetic term (HD with a = 4.071, uj = 2.4575, N = 400, L = 27r and 
h = l/407r. Both are averaged over 100 steps around each time step. The 
absorbing boundary is set at p = zizSn, and the ballistic motion occurs in the 
blue region in the schematic phase space depicted in Fig. [2]||a). (b) Long time 
behavior of P^°'^"^ for various L. 



The parameter L controls the border from which the ballistic propagation begins to 
occur. Furthermore, to avoid the recurrence of the wavepacket to the initial domain, 
we put the absorbing boundary at \p\ = PcntoS- More precisely, a projection operator 
P, which satisfies {p\P\i^) = for \p\ > Pcutos, and {p\P\i') = {p\i)) for \p\ < Pcutos, is 
applied in each time step. After penetrating through the torus region, the wavepacket 
propagates in chaotic seas up to \p\ ~ L. Then it moves ballistically in the region \p\ > L 
and is absorbed at \p\ = PcutoS- Note that the original map (2) is recovered in the limits 
L — >• oo and setting di = d2 = 0. 

We launch the wavepacket from the center of the torus region po = and observe 
the same quantity P^°^'^^ introduced above. As shown in Fig. [3l the enhancement of 
tunneling amplitude is again remarkable as compared to the system without outside 
ballistic propagation. We further notice that, as seen in Fig. [3](a), although the 
amplitude of wavefunction on the torus region is gradually reduced as a result of 
tunneling leakage, it is still distinguishably localized on the torus region. Figure [3](b) 
also exhibits that p*°''"'^ remain finite after exceedingly long time steps. Especially in 
cases of L = 27r and L = An, even a signature of saturation can be detected. 

To see that the observed process is certainly a chaos-involved one, in other words, 
not only the localization length but also the nature of chaos, especially the strength of 
chaos controls the transition amplitude, we change the nonlinear parameter K in our 
system. Taking into account that the localization length is proportional to [21], we 
plot in FigHlthe tunneling rate 7 as a function of K for several scaled values L = L/K"^. 
In this experiment, the initial wavepacket is again placed inside the torus region. The 
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Figure 4. The tunneling rate 7 as a function of the kicking strength K. 7 is 
evaluated by fitting the initial decay of P^^'^'^ for n < 10^ or 10^ depending on 
K. The length L of the chaotic region without noise (see Figl2^a)) is scaled as 
L = L/K^. Here Pd = O.Gtt, a = 4 + O.lVs and to = 27r(l + V7)/W. 

kinetic term is replaced by 

1 1 

T'{p) = -a{p - Pd) + uj + -a{p - Pd) tanh b{p - Pd) , (9) 

which gives phase space whose lower part {p < 0) is all covered with KAM circles. The 
kinetic term ([9]) is obtained by letting the second pd ^ 00 and putting di = —apd + oj. 
Such manipulation is necessary to avoid the wavepacket penetrating through the KAM 
domains and reaching the opposite chaotic region. The result shows that the degree of 
enhancement depends on the strength of chaos. Our interpretation for the result is that 
the density of stable and unstable manifolds becomes large with the increase of K, the 
number of complex paths responsible for the tunneling transition increases, which yields 
large tunneling amplitude. 

These numerical tests, together with semiclassical analysis using complex 
trajectories show that chaotic tunneling and dynamical localization are strongly 
correlated. If the interference yielding dynamical localization is destroyed, potentially 
existing tunneling trajectories show up, which leads to the drastic enhancement of 
tunneling amplitude. This makes a sharp contrast to the system coupled with the 
heat bath, in which quantum tunneling is suppressed [22]. We can alternatively say 
that torus states are sustained by surrounding dynamical localization, otherwise they 
cannot stay localized on the torus. Hence, amphibious states found in [14, 15] can be 
interpreted exactly as the flooding of chaotic tunnehng trajectories. 

An important nontrivial question still not clarified is the origin of life time of torus 
states or the tunneling rate between regular and chaotic regions. It is not obvious what 
dynamical information is needed to specify it. In the semiclassical argument, Lyapunov 
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exponents and the topological entropy of outside chaotic regions again become necessary 
ingredients since, as mentioned in the first part, dominant tunneling orbits are controlled 
by unstable periodic orbits in chaotic regions. 

It is natural to assume that inner torus states are more robust than the outer ones 
and the life time of the former is much longer than the latter [15]. This would be 
a qualitative explanation for slowing down phenomenon, but the problem looks more 
subtle. As shown in Fig{Tt^c), the orbit shown by the red curve is launched at an inner 
torus as compared to the orbits shown by green and yellow curves (see ^ coordinate in 
FiglU^a) and also note that the initial condition is place on pq = 0). Nevertheless, it 
gain smaller imaginary action ImSn, since it approaches the real plane directly, whereas 
the latter two orbits launched at outer tori take side trips, which cause additional gains 
of ImSn- This suggests the inner torus has a smaller life time than the outer one and 
so the tunneling rate does not necessarily follow a simple order even in a clean setting 
as the present model. 

It is also necessary to examine how the life time of (complex) classical orbits inside 
the torus region is related to the tunneling rate. The initial manifold, {{po = 0,C,+iTi) G 
C^}, represents the support of semiclassical wavefunction. It is found that the manifold 
starting at the region where the complex KAM domain dominates has a very long life 
time which has an entirely complex classical origin [24]. 
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